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ABSTRACT 

We use the S AURON and GMOS integral field spectrographs to observe the active galactic 
nucleus (AGN) powered outflow in NGC 1266. This unusual galaxy is relatively nearby (D=30 
Mpc), allowing us to investigate the process of AGN feedback in action. We present maps of 
the kinematics and line strengths of the ionised gas emission lines Ha, H/3, [OIII], [OI], [Nil] 
and [SII], and report on the detection of Sodium D absorption. We use these tracers to explore 
the structure of the source, derive the ionised and atomic gas kinematics and investigate the 
gas excitation and physical conditions. NGC 1266 contains two ionised gas components along 
most lines of sight, tracing the ongoing outflow and a component closer to the galaxy systemic, 
the origin of which is unclear. This gas appears to be disturbed by a nascent AGN jet. We 
confirm that the outflow in NGC 1266 is truly multiphase, containing radio plasma, atomic, 
molecular and ionised gas and X-ray emitting plasma. The outflow has velocities up to ±900 
km s _1 away from the systemic velocity, and is very likely to be removing significant amounts 
of cold gas from the galaxy. The LINER-like line-emission in NGC 1266 is extended, and 
likely arises from fast shocks caused by the interaction of the radio jet with the ISM. These 
shocks have velocities of up to 800 km s _1 , which match well with the observed velocity of 
the outflow. Sodium D equivalent width profiles are used to set constraints on the size and 
orientation of the outflow. The ionised gas morphology correlates with the nascent radio jets 
observed in 1 .4 Ghz and 5 Ghz continuum emission, supporting the suggestion that an AGN 
jet is providing the energy required to drive the outflow. 

Key words: galaxies: individual: NGC 1266 - ISM: jets and outflows - galaxies: jets - galax- 
ies: elliptical and lenticular, cD - galaxies: evolution - galaxies: ISM 



1 Introduction 

In recent years the idea feedback from an active galac- 
tic nucleus (AGN ; e.g. ISpringel. Pi Matteo & Hernquistl 120051 : 
ICroton et al]|2006h could be responsible for the quenching of star 
formation has grown in popularity. Such quenching seems to be 
requ ired to create the red-sequence galaxies we observe today 
(e.g. iBaldrv et alj|2004l) . There is circumstantial evidence to sup- 
port A GN-driven quenching, such as the study bv lSchawinski et al] 
suggesting that AGN are predominantly found in green val- 
ley galaxies, but direct evidence for removal/heating of cold star- 
forming gas is rare. 

The physical mechanism by which an AGN could drive 
molecular gas out of a galaxy is still debated. Radiation pres- 
sure is thought to be important in star fo rmation-driven outflows 
(e.g. iMurrav. Ouataert & Thompson! l200l) . and is pote ntially im- 
plicated in AGN-powered 'quas ar mode' outflows (e.g. lArav et alj 
Il999l ; iKurosawa & Progj|2009h . Kinetic feedback from an AGN 
jet can provide sufficient power to directly p ush through the ISM 
of a galaxy and entrain o r destroy it (e.g. iRosario et alj l20ld : 
iMcNamara & Nulserj|2012h . but it is unclear if the geometry of 
a bipolar jet, which often emerges perpendicular to the nuclear 
disk, will allow the jet to remove the ISM from an entire galac- 
tic disk. Broad-line winds can deposit significant momentum into 
gas surrou nding an AGN, wh ich could also lead to large out- 
flows (e.g. lOstriker et ai1l20ld) . Alternatively, heating by X-rays 
and cosmic rays could destroy/alter the molecular clouds close to 
an AGN, removing the need to expel them from the galaxy (e.g. 
iBegelman. de Kool & Sikordl99ll : lFerland et alj2009l> . These pro- 
cesses should be distinguishable if we can identify and study local 
galaxies where AGN feedback is ongoing. 
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Our recent Combined Array for Research in Millimetre-wave 
Astronomy (CARMA) and Sub-Millimetre Array (SMA) obser- 
vations of the nearby lenticular galaxy NGC 1266 suggest that 
it harbors a massive AGN-driven molecular outflow, providing 
an excellent local l aboratory for studying AGN-driven quenching 
lAlatalo etalJl201lL hearafter A2011). NGC 1266 is a nearby (D= 
29.9 Mpc; derived from recession velocity in lCappellari et al 1201 ll : 
hereafter ATLAS 3D Paper I), early-type galaxy (ETG) in the south- 
ern sky (8 = —2°), which was studied as part o f the ATLAS 3D 
proje ct. A three colour image of this galaxy (from lKennicutt et alj 
120031) is presented in Figure [JJ While typical CO spectra from 
early-type galaxies reveal the double-horned profile characteris- 
tic of gas in a relaxed disk with a flat rotation curve, the spec- 
trum of NGC 1266 shows a narrow central peak (FWHM ~120 
km s _1 ) with non-Gaussian wings o ut to at least ±40 km s _1 
with respect to the systemic velocity dYoung et alj|201 lb . Imaging 
of the high-velocity components using the SMA revealed that the 
wings resolve into redshif ted and blueshifted lob es (A2011), coin- 
ciden t with Ha emissio n l Kennicutt et alj[2003l) . 1.4 GHz contin- 
uum teaan & Kl6cknerll2006h 7 and thermal bremsstrahlung emis- 
sion (detected with Chandra; A2011; Fig. 3). Molecular gas ob- 
servations suggest that 3xl0 8 Mq of molecular gas is contained 
within the central 100 pc of NGC 1266, and that at least 5x 10 7 
Mq of this gas is involved in a molecular outflow (A201 1). This is 
thus the first observed large-scale expulsion of molecular gas from 
a non-starbursting ETG in the local universe, and this presents a 
unique opportunity to study this powerful process in action. 

In this paper we present SAURON (Spectrographic Areal 
Unit for Research on Optical Nebulae) and Gemini Multi-Object 
Spectrograph (GMOS) integral-field unit (IFU) observations of the 
ionised gas in NGC 1266. By investigating the ionised gas kine- 
matics and line ratios we hope to constrain the outflow parameters 
and ionisation mechanisms and thus shed light on the mechanism 
driving gas from the galaxy. In Section[2]we present the data, and 
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Figure 1. SINGs jKennicutt et alj2003l) B,V and R band composite three 
colour image of SO galaxy NGC 1266. The white bar shows a linear scale 
of 1 Kpc (6^94 at an adopted distance of 29.9 Mpc; ATLAS 3D Paper I). 
Overlaid are the total field of view of our SAORON IFU (red) and GMOS 
IFU (blue) observations. 



describe our reduction processes. We then present the derived maps 
of the gas kinematics and line fluxes. In Section [3] we discuss the 
kinematic structure of the system, gas excitation mechanisms and 
the driving force behind the outflow. Finally we conclude and dis- 
cuss prospects for the future in Section[4] 

2 Data Reduction and Results 
2.1 SAURON data 

SAURON is an integral-field spectrograph built at Lyon Observa- 
tory and mounted at the Cassegrain focus of the Willi am Herschel 
Teles cope (WHT). It is based on the TIGER concept dBacon et alj 
1 19951) . using a microlens array to samp le the field of view . De- 
tails of the instrument can be found in iBacon et alj feOOlh . The 
SAURON data of NGC 1266 was taken at the William Herschel 
Telescope (WHT), on the night of 10-11 January 2008, as part of 
the ATLAS 3D observing campaign (ATLAS 3D Paper I). The galaxy 
was observed with the low-resolution mode of SAURON, covering 
a field of view of about 33" x 41" with 0"94 x 0"94 lenslets. The 
field of view (FOV) of our observations is shown in red on Figure 
Q] SAURON covers the wavelength range from 4810-5350 A with a 
spectral resolution of 105 km s _1 . 

The basic reduction of the SAURON observation was accom- 
plished using the standard ATLAS 3D pipeline. Details of this pro- 
cess, including extracti on of the stellar kinem atics are presented 
in ATLAS 3D Paper II dKrainovic eTall 1201 lh . In brief, the two 
observed datacubes w ere merged and processed as described in 



lEmsellem et alj 1 20041). using the V oronoi binning scheme devel- 
oped by Cappellari & Copinf d2003l) . This binning scheme max- 



imises the scientific potential of the data by ensuring a minimum 
signal-to-noise ratio of 40 per spatial and spectral pixel. This does 
however result in an non-uniform spatial resolution, here varying 
from 0'.'8 x 0'.'8 for unbinned spaxels in the central regions, to 
10" x 7"in the largest outer bin. 

The SAURON stellar kinematics were derived using a penal- 



ized pixel fitting routine dCappellari & Emsellem|[2004l) . providing 
parametric estimates of the line-of-sight velocity distribution for 
each bin. Du ring the extraction of the stellar kinematics, the GAN- 
DALF code dSar zi et al. 20o3) was use d to simultaneously extract 
the ionised gas line fluxes and kinematics. The standard GANDALF 
reduction completed in the pipeline (using a single gaussian for the 
lines) is insufficient in this source, due to the complex structure of 
the ionised gas outflow (see Figure(2]l. We have reanalyzed the dat- 
acube using a multi-gaussian technique (as described below) after 
the subtraction of the stellar continuum. 

2.1.1 Emission line fitting 

The SAURON spectra include the H/3, [OIII] and [NI] ionised gas 
spectral lines. As can be seen in Figure[2]some of the binned spax- 
els show clear signs of having two ionised gas components along 
the line of sight with different velocities. In order to fit these pro- 
files we created an IDL pr ocedure based on the Non-Linear Least 
Squares Fitting code mpfit dMarkw ardt 2009). In this procedure we 
perform two fits, and compare the chi-square to determine if two 
components are needed at each position. 

In the first fit, we assume a single ionised gas component is 
present, and fit the H/3, [OIII] and [NI] lines with single gaussians. 
These gaussians are constrained to have the same kinematics (ve- 
locity and velocity dispersion). Additionally, we confine the veloc- 
ity of the lines to be within 1000 km s _1 of the galaxy systemic 
(2170 km s _1 ; ATLAS 3D Paper I), and to have a velocity dispersion 
greater than the instrumental resolution, and less than a convolved 
velocity dispersion of ~200 km s _1 . Initial guesses at the ionised 
gas kinematics were made by assuming the ionised gas co-rotates 
with the stars, with a velocity dispersion of 120 km s . We con- 
strain the fitting by forcing each gaussian to have a peak at least 3 
times larger than the noise in the continuum, or to be zero. Any flux 
which had a la error bar that included zero was set to zero. Initial 
guesses of the line fluxes were estimated by taking the maximum 
flux within the allowed velocity range of each line. The two [OIII] 
lines in our spectral range have a fixed line ratio determined by the 
energy structure of the atom, and we fixed the line ratio assumed 
in our fit to agree with the observed line ratio (F 5qq7=2.99xF4959: 
IStorev & Zeipperj|200ri iDimitriievic et alj|2007t) . 

In the second fit we assume two, independent ionised gas com- 
ponents are present in each bin, and fit each component with its own 
set of independent linked gaussians. As before, it is assumed that the 
lines in each component trace the same kinematics (velocity and ve- 
locity dispersion). Once again we confine the velocity of each of the 
components to be within 1000 km s _1 of the galaxy systemic, and 
to have a velocity dispersion greater than the instrumental resolu- 
tion. We use different upper bounds for the first and second com- 
ponents of the gas distribution. Component one is forced to have a 
velocity dispersion less than 200 km s _1 , as before. In general the 
second component is needed where the outflow is present, and is 
thus allowed to have a higher velocity dispersion. We allowed the 
lines in the second component to have a maximum velocity disper- 
sion of 360 km s _1 . In practice however good fits were found with 
velocity dispersions <300km s _1 . The same limits were used on 
the line fluxes as described above. 

Once the two fits described above were complete for each bin, 
we tested (using an F-test, as implemented in the mpfit package 
Markwardt 2009) if adding the additional free parameters to our 
model of the emission lines produced a significantly better fit, over 
and above the improvement expected when one adds free param- 
eters. The F-test can be used as an indicator of where fitting two 
components produces better models, but the best threshold to take 
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Figure 2. Stellar emission subtracted SAURON spectrum from single bins 
(black solid line) with line identifications. The top panel shows the spaxel 
with the largest line flux (x=0", y=-4"), the middle panel shows the spaxel 
with the biggest difference in the fitted velocities (x=0"8 y=-5'.'17) and the 
bottom panel shows the lowest flux region where a two component fit can be 
constrained (x=-4'.'0 y=-5" 17). Overlaid is the two component fit produced 
by out fitting routine, as described in Section l2. 1.11 Component one (nearest 
the galaxy systemic) is shown in orange, while component two is shown in 
blue. The red dashed line shows the sum of component one and two, which 
closely matches the observed data. 



should be determined by visually inspecting the fits obtained (as 
the values tested for are at the extreme edge of the possible distri- 
bution). In this work we visually chose a threshold that corresponds 
to an improvement in the chi-square of 60% when adding in the 
additional parameters. When a spaxel did not satisfy this criterion 
then the values from the single gaussian fit were used, and the sec- 
ond component set to zero. Where two components were found to 
be necessary we denoted the component closest to the galaxy sys- 
temic as component one, and the faster component as component 
two. 

In an attempt to ensure that the fits were robust, and spatially 
continuous, we implemented an iterative fitting regime where the 
fitting processes described above were performed for each spaxel 
in turn. Then the resulting two dimensional flux and velocity maps 
were smoothed using a gaussian kernel, and then these smoothed 
values used as the initial guesses for the next iteration of the fit- 
ting procedure. Using this procedure we found that the parameters 
usually converged within three iterations, with very little variance 
between fitting attempts. Figure [2] shows SAURON spectra from 
a single bin, overlaid with the two component fit. The top panel 
shows the spaxel with the largest line flux, the middle panel shows 
the spaxel with the biggest difference in the fitted velocities and the 
bottom panel shows the lowest flux region where a two component 
fit can be constrained. Clearly in the low flux regions of the cube the 
fitted velocities are driven by the [OIII]5oo7 line, and the parameters 
have correspondingly higher uncertainties. 

In Figure [3] we show the observed kinematics for component 
one (panels a & b) and component two (panels c & d). The veloc- 
ity dispersion maps have had the instrumental dispersion quadrati- 
cally subtracted. We also show the stellar velocity field derived from 
these SAURON observations (panel e) for reference, as presented in 
Paper II . In Figure|4]we show the fitted line fluxes. 

2.2 Gemini GMOS data 

In addition to the SAURON data, we obtained Gemini GMOS-IFU 
observations of the central parts of NGC 1266, providing higher 
spatial resolution and a longer wavelength coverage. The GMOS 
IFU uses a lenslet array of 1500 ele ments to feed individual po- 
sitions on the sk y to optical fibres dAllington-Smith et al]|20ol 
iHook et al.l2004l) . The total field of view of the IFU is 5" x 7", with 
a spatial sampling of df.'2. The Gemini GMOS-IFU observations 
of NGC 1266 were taken over the nights of 24th, 26th and 27th of 
January 2009 at the Gemini North telescope (program GN-2008B- 
DD-1). We used a four point dither pattern to extend our coverage 
to a total field of view of ~9'.'l x 12'.'5, around the optical centre 
of the galaxy. The resultant field of view (FOV) of our observa- 
tions is shown in blue on Figure [T] The low resolution R150 grat- 
ing was used, resulting in a spectral resolution of ~185 km s _1 (at 
6500 A) over the wavelength range 5000 - 7300 A. Two different 
blaze wavelengths (688 and 700 nm) were used on different expo- 
sures to allow continuous spectral coverage by averaging over chip 
gaps/bridges. 

In order to reduce the GMOS IFU data we uti lized a data re- 
duction pipeline, as used in Ivan de Ven et alj ll2010h . This pipeline 
calibrates and flat fields the data, before it is trimmed and resam- 
pled into a homogeneous data cube. This cube was binned using the 
Voronoi binning technique o f ICappellari & Copinl ( Eoolh . ensuring 
a signal-to-noise ratio of 40 per spatial and spectral pixel. Due to 
the low spectral resolution and the depth of the exposures we de- 
tect line emission to high significance over almost the entire IFU 
cube, but were unable to detect stellar absorption features to high 
significance. As we are interested in the ionised gas kinematics in 
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Figure 3. Ionised gas kinematics derived from the SAURON IFU data reduction process discussed in Section l2.1.1l In the top row (panels a and b) we display 
the kinematics of component one (confined to be closest to the galaxy systemic velocity). Bins where only one ionised gas component is required are also shown 
in component one. The kinematics of the faster component is shown in the middle row (panels c and d). The ionised gas velocity is displayed in the left panels 
(a,c), and the velocity dispersion in the right panels (b,d). The stellar velocity field of this galaxy derived from these same observation in ATLAS JD Paper II is 
shown as a comparison in panel e, with stellar continuum flux contours overlaid. These plots are centred around the galaxy position given in ATLAS 3D Paper 
I. 



this work we simply wish to remove the (relatively smooth) stel- 
lar continuum. We do this by fit ting the stellar continuum with the 
penalized pixel fitting routine of Capp ellari & Emselleml d2004l) . as 
used for our SAURON data. We were able to constrain the number 
of stellar templates required using our best fit to the SAURON data. 
Once the stellar continuum was successfully removed we were left 
with a cube containing the ionised gas emission only. 

The spectral range of our cube includes various ionised gas 
emission lines. The [OIII] and [NI] are included in our GMOS spec- 
trum in the region which overlaps with the SAURON spectral range. 
These lines however appear very weak because they are at the edge 
of our band pass, where the throughput is low. The main strong 
lines we detect are Ha, [NII]6548,6583 and [OI]e3oo, while Hel and 
[Nil] 5754 are detected more weakly. We choose to fit the kinemat- 
ics of the gas emission on the strong lines only, and impose these 



kinematics when measuring the fluxes of weaker lines. Example 
specta extracted from the cube are shown in Figure [5] We show 
here the region around the Ha, [Nil] and [SII] lines only. These 
specta were selected to lie at the spatial position with the highest 
line flux (top panel), the spaxel with the largest difference between 
the two fitted velocities (middle panel) and at the lowest flux bin 
in which we can constrain two components (bottom panel). With 
the low spectral resolution of this data the lines are blended, how- 
ever we clearly require two components to fit the line emission. We 
also detect Sodium D (NaD) absorption against the stellar contin- 
uum (an example spectrum is shown in Figure |6). We describe the 
fitting procedure used for the gas emission lines in detail in Section 
12.2.11 and the procedure used to measure the parameters of the NaD 
absorption in Section l2.2.2l 
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Figure 4. Ionised gas line fluxes derived from the SAURON IFU data reduction process discussed in Section 12.1.11 In the left column we display the flux of 
component one (confined to be closest to the galaxy systemic). Bins where only one ionised gas component is required are also shown in component one. The 
flux of the faster out-flowing component is shown in the right hand column. The top row shows the H/3 line fluxes (panels a & b), the second row the [OIII] 
line fluxes (panels c & d), and the third row the [NI] fluxes (panels e & f). Fluxes are in units of 10 — 16 erg s — 1 cm -2 arcsec - 2 in each SAURON bin. 



2.2.1 Emission line fitting 



GMOS emission line fitting was carried out as described in Section 
12.1.11 with some modifications, described below. We fit here the 
Ha, [Nil], [SII] and [01] lines with single and double gaussians. 
Initial guesses for the velocity and velocity dispersion were made 
using the derived kinematics from the SAURON cube. The two 
[Nil] lines in the Ha region of the spectrum have a fixed line ratio 
determined by the energy structure of the atom, and we forced the 
line ratios to the theoretical line ratio (F6584 =2.95xF6548: Acker 
1989). The [SII] doublet is an electron density tracer, with the 
line ratio F6731/F6717 varying from .459 in the high density limit 
to 1.4 3 at the low density limit (e.g. |Pe Robertis. Dufour & Hurra 
Here we constrain the [SII] lines ratio to lie somewhere 
within this region. We here allowed the lines in the second com- 
ponent to have a maximum velocity dispersion of 500km s _1 , but 



again note that good fits were found in most bins with velocity dis- 
persions < 300km s _1 . 

As for the SAURON data, to ensure that the fits were robust, 
and spatially continuous we implemented an iterative fitting regime 
where the fitting processes described above were performed multi- 
ple times, using a smoothed version of the output from the previous 
run as the initial conditions. Here we found that the parameters usu- 
ally converged within four iterations, again with very little variance 
between fitting attempts. Figure|5]shows the GMOS spectra from a 
single bin (from the same spatial region as selected before), overlaid 
with the two component fit, as found by this procedure. 

In Figure [7] we show the observed kinematics for component 
one (panels a & b) and component two (panels c & d). The velocity 
dispersion maps have had the instrumental dispersion quadratically 
subtracted. In Figure[8]we show the fitted line equivalent widths for 
the strong lines, with fitted component one in the top row and com- 
ponent two in the bottom row. We calculate the equivalent width in 
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Figure 5. Stellar subtracted GMOS spectrum (black solid lines)from single 
bins with line identifications. Shown in the top panel is the bin with the 
spaxel with the peak line flux (x=0'/91, y=-2'.'08), the middle panel is the 
spaxel with the largest difference between the two fitted velocities (x=-0 / . / 2, 
y=-4'.' 29), and the bottom panel shows the bin with the lowest total flux in 
which we are able to fit two components (x=-l'.'85, y=-5"73). Overlaid is the 
two component fit produced by out fitting routine, as described in Section 
12.21 Component one (nearest the galaxy systemic) is shown in orange, while 
the faster component two is shown in blue. The red dashed line shows the 
sum of component one and two, which closely matches the observed data. 
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Figure 6. GMOS spectrum of the NaD region from a single bin (x=-l"0, 
y=-l'.'8; black solid line). Overlaid is the two component fit produced by our 
fitting routine for the emission lines, as described in Section l2~2l Component 
one (nearest the galaxy systemic) is shown in orange, while the out-flowing 
component two is shown in blue. The red line shows the sum of component 
one and two (and our fit to the Nal sky line), which closely matches the 
observed data. Shown in green is our fit to the NaD absorption trough. 

the standard way by finding the width of a rectangle, with a height 
which is the same as the average stellar continuum flux in the region 
of the lines, which has the same area as the observed lines. As the 
GMOS data were taken in non-photometric conditions, we will use 
only ratios of the line fluxes from this point on. 

2.2.2 Absorption line fitting 

The sodium absorption lines at 5890 A and 5896 A are detected 
in our GMOS IFU data (Figure [6](. This feature is unlikely to be 
due to an imperfect stellar template leaving negative residuals after 
subtraction from the GMOS spectrum, as the fitted velocities of the 
absorption feature do not match the stellar velocities assumed when 
fitting the template. 

In order to extract the absorption depths, and determine the 
neutral gas kinematics we jointly fit the absorption doublet, and the 
neighbouring Hel emission line (and Nal skyline). We fix the ve- 
locity and velocity dispersion of the Hel line using the best solution 
for each bin derived from the stronger lines (as described in Sec- 
tion [2T27TJ- We then fit this line with two gaussian components, as 
before, to determine the line fluxes. 

Simultaneously we fitted the NaD absorption doublet (which is 
blended in our data), assuming a gaussian profile for both lines. For- 
mally an absorption line should be fitted with a Voigt profile, as the 
absorption has an intrinsic Lorentzian shape, which has been con- 
volved with the instrumental gaussian response. In the low spectral 
resolution data we present here however we fit gaussian profiles, as 
the instrument response function is much broader than the intrinsic 
absorption. If one fits a Voigt profile to our data, the best fit profiles 
always tend towards a pure gaussian, with a Lorentzian width (r) 
of zero, validating such an approach. We do not fix the NaD veloc- 
ity and velocity dispersion, as the absorption arises from a different 
gas phase, which may have different kinematics (see Section l3?2t . 
The velocities of the NaD hosting gas are constrained in our fit to 
lie within 1000 km s _1 of the systemic velocity, and the velocity 
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dispersion of this component is constrained to be greater than the 
instrumental, but less than 500km s _1 . At the spectral resolution 
of our data we only require a single neutral gas component in all 
spaxels in order to fit the absorption profiles well. 

In Figure[9]we show the observed absorption equivalent width, 
and the kinematics for the NaD absorbing gas. We calculate the 
equivalent width in the standard way, as above. 

3 Discussion 

3.1 Ionised gas emission line kinematics 

Figures [3] and UJ show the ionised gas kinematics derived from our 
multi-gaussian fitting procedure. The SAURON data has a much 
wider field of view, providing insight into the large scale kinematics 
of this object. The GMOS IFU data zooms in to the central portions 
of this object to reveal the inner regions. 

Panel a of Figure [3] shows the SAURON ionised gas kinemat- 
ics for the component nearest the galaxy systemic velocity (which 
we will hereafter call the systemic component). This component in- 
cludes gas out to a radius of ~10" (1.5 Kpc). A201 1 discussed the 
molecular gas distribution, and find evidence for a rotating molecu- 
lar disk, as well as a molecular outflow. The origin of this systemic 
component, and its relation to the molecular disk is unclear. This 
component could be due to unrelated gas components at different 
locations along the line of sight, or it may be a coherent structure 
which has been disturbed. 

Some degree of symmetry appears to exist in the gas distribu- 
tion around a line inclined ~20° from East. This may suggest some 
bulk rotation, with a kinematic position angle of «30°. If this com- 
ponent is rotating, then comparison with the stellar rotation (shown 
in the bottom row of Figure[3} shows that it is misaligned from the 
stellar rotation by ~90°, suggesting it could be in the polar plane. 
In the inner ~6"however the velocity field changes sign, and is dis- 
turbed. 

Figure |7] shows the GMOS zoomed in view of the centre of 
NGC 1266. The same disturbed features present in the SAURON 
data are observed in the GMOS velocity field (Panel a of Figure 
UJ. These features persist whatever our choice of initial velocities 
for the fitting procedure, and are located at the same spatial loca- 
tion as the most blue- and red-shifted gas in component two. When 
fitting multiple components to observed velocity profiles its always 
possible that such reversed sign velocity structures are a result of 
a mis-fitting or over-fitting of the line components. Here however 
we find the same structure from both the SAURON and GMOS data 
independently. Our iterative fitting procedure described in Sections 
12. 1 . 1 I and l2.2. H is designed to avoid discontinuous fits, and thus tries 
to remove such disturbed structures, but is unable to find better fits 
in these spaxels. The middle panel of Figures [2] and [5] show the 
emission line spectrum in the bin with the largest blue-shifted ve- 
locity in each dataset (where the disturbed systemic component is 
also detected). Clearly the outer edges of each blended line or line- 
complex (which drive our determination of the relative velocities of 
the two components) are well fit, increasing our confidence that this 
structure is real. 

If the systemic component is a misaligned rotating structure 
then we speculate that the disturbed structure visible in the inner 
parts could be where the molecular outflow detected in A2011 has 
punched through, and material is flowing inward to fill the vacuum. 
Given the dynamical timescale of gas at this radius (~30-40 Myr; 
A2011) it is difficult to imagine that the disturbed gas in the cen- 
tre of this object could exist for long without being smeared out. 
This suggests that either this feature is very young, or we are in- 



deed observing unrelated clouds along the line of sight, which are 
not dynamically linked. A2011 estimate the age of the molecular 
outflow in this object as ~2.6 Myr, so a recent cause of this fea- 
ture is not completely ruled out. Alternatively if these are unrelated 
clouds long the line of sight, the blue features south of the nucleus 
may be directly related to the outflow. Close correlation between 
some of these features and the radio jet support this conclusion (see 
Section [3~5l >. Higher spatial resolution observations would allow us 
to disentangle these two possibilities. 

Panel c of Figure [3] shows the global SAURON view of the 
gas kinematics where a second component was required. This gas 
appears to be in a two lobed structure, orientated approximately 
North- South (misaligned fro m the kinematic axis of the stars by 
«70° dKrainovic et alj|201lh ). The gas in this component has ve- 
locities up to ~800 km s _1 away from the systemic velocity. We 
denote this component as the outflow from this point. Panel c of 
Figure [7j shows this component in more detail. With the finer spa- 
tial sampling of the GMOS-IFU map we are able to detect gas with 
velocities up to ~900km s~ away from the systemic. As argued 
in A2011 the escape velocity in the centre of this object is at most 
~340km s _1 supporting the idea that the outflow in this system 
allows material to escape the galaxy. The molecular gas in the out- 
flow of this object reaches velocities of up to ~480 km s _1 , with 
a slightly smaller spatial extent (see Figure [Tot, This could suggest 
molecular gas is destroyed as it flows out of the galaxy, or that our 
observations were not sensitive to detect emission from the molec- 
ular gas at the highest velocities. 

The molecular gas in NGC 1266 dominates the total mass of 
the ISM (with a mass of ~2x 10 9 Mq), and is contained within the 
inner most ~100 pc of the galaxy (A201 1). It is very hard to explain 
how this gas lost its angular momentum if it was already present 
within the galaxy. Deep optical imaging shows some minor signa- 
tures that could be due to recent disturbances, but no signatures of 
a stellar bar (Due et al., in prep). A minor merger could explain the 
compactness of the gas, if the merger happened with the correct ini- 
tial parameters to leave the gas with little or no angular momentum. 
The dynamical mass (M1/2) of NGC1266 within a sphere of radius 
r i/2 (containing half of the galaxy light) is 1 x 10 Mq (Cappellari 
et al., in prep). A minor merger with a smaller galaxy could provide 
the gas we see, and with a stellar mass ratio of ~5:1 (assuming the 
smaller galaxy was gas rich, with a gas fraction of one) may not 
leave visible traces in the luminosity weighted galaxy kinematics. 

In fbavis et al U20H we analyzed the kinematic misalignment 
of the ionised g as in ATLAS 3D galaxies (extending the work of 
ISarzi et al]|200m in order to gain clues about the origin of the gas. 
In that work we suggested that gas with kinematic misalignments 
>30° almost certainly have externally accreted gas. If the systemic 
component is rotating (and its misalignment from the stellar ro- 
tation axis is not caused by orbit bran ching or similar; fpfennigeil 
ll985l : IContopoulos & Magnenad[l985h then this would once again 
argue for a recent minor merger/accretion. 

3.2 NaD absorption 

The sodium absorption doublet at 5890 A and 5896 A provides a 
good probe of the cold ISM in the outflow. The ionisation potential 
of Nal is lower than that of hydrogen, at only 5.1 eV. The pho- 
tons that ionize Nal are thus in the near- UV (A w 2420 A). These 
lines primarily probe the ISM in the warm atomic and cold molecu- 
lar phases, where t here is sufficien t dust extinction to allow neutral 
sodium to survive (Spit ze3ll978l) . For NaD lines to be observed, 
only relatively modest optical depths and HI column densities are 
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Figure 7. Ionised gas kinematics derived from the GMOS IFU data reduction process discussed in Section 12.21 In the top row (panels a & b)we display 
the kinematics of component one (confined to be closest to the galaxy systemic). Bins where only one ionised gas component is required are also shown in 
component one. The kinematics of the faster out-flowing component are shown in the bottom row (panels c & d). The ionised gas velocity is displayed in the 
left panels, and the velocity dispersion in the right hand panels. 



required, which makes these lines a sensitive probe of the outflow- 
ing (or inflowing) neutral ISM. 

As shown in Figure [9] we detect NaD only in the central re- 
gions, and southern part of the galaxy. As absorption lines are 
viewed against the stellar continuum, if blue shifted velocities are 
observed it is clear that the gas is outflowing, rather than inflow- 
ing. The gas kinematics show that it is likely caught in the outflow, 
with typical blue shifted velocities of ~-250km s" 1 and extreme 
velocities detected out to —500 km s _1 (which is well beyond the 
escape velocity). The position angle of the outflow traced in NaD 
(and molecular gas; see Figure [TUt is slightly different than that 
traced by the ionised gas, by ~35±5° . It is unclear if the difference 
between the two gas phases is significant. We discuss this further in 
Sectionl331 

The presence of NaD in the outflow provides further evidence 
that NGC 1266 hosts a multi-phase outflow, and that cold gas is 
being removed from the galaxy. As discussed above, the outflow 
extends to higher velocities when traced by ionised gas than when 
using a dense gas tracer. Either the gas is dissociated further out in 
the outflow, or we no longer have the sensitivity to detect it. 

To observe the NaD lines, the extinction must be sufficient for 
the optical depth (t) to be > 1 at 24 20 A which corresponds to 
an A v <; 0.43 mag in the V-band (for a lCardelli. Clayton & Mathisl 



Il989l extinction law) and to a HI column density of >8xl0 
cm~ 2 . The NaD lines we observe are likely to be associated with 
the HI in this source, which is detected in absorption by A2011. 
They find a total HI column density of Nhi = 2.1 x 10 21 cm~ 2 in 
front of the continuum source in NGC 1266, and estimate that the 
HI column depth in the outflow is ~ 8.9xl0 20 cm -2 , consistent 
with our detection of NaD. 

Our observations provide an alternative way to estimate the 
reddening, and thus the atomic gas column density in each spaxel. 
The equivalent width (EW) of the NaD absorption lines has been 
shown to correlate with reddening. Here we use the relation of 
iTuratto. Benetti & Cappellarol d2003h . derived from low resolution 
spectroscopic observations of supernovae: 

E ( B zlQ. = _ ,oi + 0.16 x EW(NaD). (1) 
mag 

where E(B-V) is the standard colour excess. This can be combined 
with the relati on between E(B-V) and HI colu mn density, such as 
that derived bv lBohlin. Savage & Drake] dl978h : 

N(HI) = E(B - V) 
cm- 2 0.2 x 10~ 21 ' 

Using these relations we find that the average HI column density in 
the outflow of NGC 1266, as derived from the NaD EW (displayed 
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Figure 8. Ionised gas line equivalent widths derived from the GMOS IFU data reduction process discussed in Section l2~2l In the first row we display the EW 
of component one (confined to be closest to the galaxy systemic). Bins where only one ionised gas component is required are also shown in component one. 
The EW of the faster out-flowing component is shown in the second row. The Ha line EWs are shown in panels a and d, the [Nil] line EWs in panels b and e, 
and the [SII] EWs in columns c and f. 




Figure 9. NaD atomic gas kinematics derived from the GMOS IFU data reduction process discussed in Section 12.2.21 In panel a we display the equivalent 
width of the absorption line. Panel b shows the derived gas kinematics, and panel c the derived velocity dispersion. 
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Figure 10. GMOS Ionised gas outflow kinematics (panel a), Hce outflow EW (panel b) and NaD kinematics (panel c), as in Figures Q [8] and [9] but plotted in 
grey scale. These are overplotted with the CO observations of A201 1. White contours show the molecular core, and blue and red contours the redshifted and 
blue shifted gas, respectively. 



in Figure[9} is (1.2±0.6)x 10 21 cm 2 , consistent with the HI ab- 
sorption spectrum presented in A201 1. 

The above estimate of N(HI) is very simplistic, and assumes 
that the NaD absorption comes from a 'screen' in front of the 
galaxy. In fact, due to the small spatial scales we are probing here, 
some of the stellar continuum detected by GMOS will come from 
stars in front of the outflow gas, decreasing the measured EW. The 
amount of contamination will change with radius. Just to the North 
of the galaxy nucleus the gas velocity is redshifted with respect to 
the galaxy systemic velocity. We are thus likely to be seeing through 
to the receding part of the outflow. If we assume that the outflow is 
symmetric, then as we only detect redshifted material at the very 
centre, this implies that the outflow is not in the plane of the galaxy 
(as otherwise the outflow on the northern side of the galaxy would 
be silhouetted by the same amount of stellar continuum). 

It is possible to use the above facts to introduce a new sim- 
ple model which can help constrain the geometry of outflows. In 
this model we treat the inner regions of NGC 1266 as a planar 
disk, with a vertical scale height of 700 pc. The exact value we 
choose for the scale height does not critically change our results, 
so here we ad opt a value at the high end of those found for the 
Milky Way bv lChen et alj {2001). As the inner surface brightness 
profile of NGC 1266 (within «5") is reasonably flat (Krajnovic et 
al., submitted) in this simple model we neglect the change in sur- 
face brightness with radius of the system. We input an outflow with 
a constant surface density, and a variable length (L), which is in- 
clined from the galaxy plane by an angle 6. From EquationQ]and[2] 
we calculate the expected EW(NaD) produced by such an outflow. 
The measured EW will however be contaminated by the luminos- 
ity of the stars between the observer and the absorbing material. 
At each radius we calculate the amount of stellar luminosity above 
and below the position of the outflow, and produce a model for the 
observed EW(NaD) profile. 

Figure [TT] demonstrates the model geometry, and includes 
three sight-lines as examples of this procedure. In Sightline (a) the 
outflowing material is in front of the majority of the stars in the 
galaxy, and thus the expected NaD EW will be relatively uncon- 
taminated. Sightline (b) passes through the galaxy centre, and 50% 
of the stars in the galaxy are in front of the absorbing material, thus 
our measured EW will be 50% smaller. Sightline (c) has almost all 
of the stars in the galaxy in front of the absorber, so the absorption 
on this sightline is likely undetectable. The exact viewing angle be- 



A A 




Figure 11. Simple geometry assumed for the NaD absorption model de- 
scribed in Section 13.21 The galaxy is represented as a flat slab with expo- 
nentially declining surface density of scale height 700 pc. The outflow has 
a constant surface density, a variable length (L), and an inclination from the 
galaxy plane of angle 9. Three sight lines (referred to in the text) are labeled 
(a,b & c). 



tween the galaxy plane and the line of sight introduces a constant 
geometric factor, which can be neglected when one normalises the 
EW profile. The profile shape does not depend on the exact col- 
umn density assumed for the outflow, and hence this factor is also 
removed when one normalises. 

Models of the above sort can be fitted to the observed NaD EW 
distribution. Figure [12] shows the observed NaD EW distribution 
(black points) extracted from our GMOS observations in a pseudo 
slit of 1" wide aligned in the North-South direction over the centre 
of the galaxy (along what seems to be the major axis of the outflow). 
Overplotted on the observed data is the best fit model for the outflow 
produced in the way described above (where additionally the model 
has been convolved with a gaussian to match the average seeing of 
our observations). Models with inclinations between and 90°, and 
with sizes between 1 and 500 pc were created, and compared to the 
data by calculating the reduced chi-squared. The statistical error is 
calculated by finding the area of parameter space around the best 
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Figure 12. The observed NaD EW distribution (black points) extracted from 
our GMOS observations in a pseudo slit of 1 " wide, aligned in the north- 
south direction over the centre of the galaxy. The solid black line is the best 
fit model for the outflow EW produced in the way described in Section l3~2l 
convolved with a gaussian to match the average seeing of our observations. 
The best model displayed here has an outflow inclination (with respect to 
the galaxy disk) of 53°, and a linear outflow size of 400 pc. The dashed line 
shows the detection limit of our observations. 



fit where the reduced chi-squared changes by 1. The best fit model 
has an inclination (between the galaxy plane and the outflow) of 
53±8°, and a size of 400±50 pc. A line corresponding to such a 
model is shown in Figure [T2l 

A201 1 estimate the outflow to have a total extent of 465 pc, 
and an inclination with respect to the plane of the sky of roughly 
20° (using the the average offset between the centroids of the out- 
flow with respect to the nucleus divided by the average extent of 
the lobes). If we use the inclination of 34±5° for the galaxy, as in 
A201 1, combined with our estimate of the outflow inclination (with 
respect to the galaxy) this leads to an estimate of the outflow incli- 
nation with respect to the plane of the sky of 18.5±10°. Both the 
inclination of the outflow and its size derived from our simple model 
are thus consistent with the estimates presented in A20 11. If we al- 
ternati vely use the inclination of 26±5° estimated in iDavis et al.l 
J20T l|), within our error bars we still get a result consistent with the 
estimate of the outflow inclination in A201 1. 

This conclusion that the outflow in NGC 1266 is orientated out 
of the plane of the galaxy is supported by the observed kinematic 
misalignment between the galaxy stellar kinematics, the ionised 
gas kinematics, and the outflow. A201 1 also argue that the outflow 
in NGC 1266 is orientated out of the plane of the galaxy, based 
on the differential reddening between the northern and southern 
lobes of ionised gas emission. In Seyfert galaxies no correlation be- 
tween the direction of AG N jets and the galactic plane exists (e.g. 
lUlvestad & Wilson! 19841) , hence the misalignment of the outflow in 
this galaxy is perhaps not unexpected. 

Two parameter models of this type employed here are highly 
simplistic, and ignore a large amount of relevant physics. Despite 
these reservations a good fit to observational data can be obtained 
(as shown in Figure [T2l l. and the estimated parameters seem reason- 
able. This highlights the power of NaD observations to constrain the 
geometry of outflows. 
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Figure 13. BPT ^Baldwin. Phillips & Terlevichlll98lh type diagrams for the 
inner parts of NGC 1 266. The Y-axis of each plot shows the [OHTJ/H/3 ratio 
derived from our SAURON data, and it is plotted versus the [NII]/H a (panel 
a; top), [SII]/Ha (panel b; middle) and [OI]/Ha (panel c; bottom) line ratios 
from our GMOS observations. In the bottom right of each plot is the typical 
error bar associated with each point, derived from the formal fitting errors 
and a monte-carlo error analysis of the fluxes returned by our fitti ng rou- 
tines. We overplot the diagnostic lines flKewlev et aI]j200lLEoo3) . which 
indicate the dominant line excitation mechanism. 
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Figure 14. Observed line ratios for spaxels in the inner parts of NGC 1266, 
plott ed over a shock m odel with a super solar metallicity (Ri 2Zq), taken 
from lAllen et alj j2008b . The Y-axis of each plot shows the [OIII]/H/3 ratio 
derived from our SAURON data, and it is plotted versus the [NII]/Ha (panel 
a; top), [SII]/Ha (panel b; middle) and [OI]/Ha (panel c; bottom) line ratios 
from our GMOS observations. In the bottom of each plot is the typical error 
bar associated with each point, derived from the formal fitting errors and a 
monte-carlo error analysis of the fluxes returned by our fitting routines. 



3.3 Gas excitation 

Figures [4] and [8] show the distribution of ionised gas in NGC 1266. 
Ionised gas is detected everywhere in the inner «10" (a projected 
distance of 1440 pc) of the galaxy, however the brightness distri- 
bution is not smooth. As seen in narrow band images (and previ- 
ously discussed in A2011) the ionised gas emission is brighter to 
the south of the nucleus, and has a distinctive 'kidney bowl' shape 
(best seen in Figure [8} • This structure correlates with the position 
of the outflow (and with the radio jet; Section [3~5l >, The asymmetric 
brightening of the lines is present in both the systemic and outflow 
components, suggesting these components are linked in some way. 
In the rest of this section we use ratios of the line fluxes to under- 
stand the ionisation mechanism powering this emission. 

Using the combination of the SAURON and GMOS 
IFU data is is possible to construct BPT diagrams 
tealdwi~ Phillips & Terlevichl Il98ll) for the inner parts of 
NGC 1266. In order to use ratios from both datasets we interpolate 
the SAURON IFU map onto the GMOS bins. As the SAURON bins 
are in general larger, some GMOS bins will have identical values 
of the SAURON line ratios. Figure [13] shows the [OIIIj/H/3 ratio 
derived from our SAURON data, plotted versus the [NII]/H Q (top 
panel), [SII]/Hq (middle panel) and [OI]/Hq (bottom panel) line 
ratios from our GMOS observations. In the bottom right of each 
plot is the typical error bar associated with each point, derived 
from the formal fitting errors and a monte-carlo error analysis of 
the different fluxes returned by our fitting routines. 

Onto these diagram s we overplot the diagnostic lines of 
iKewlev et all j200lLl2006h . which indicate the dominant line exci- 
tation mechanism. In all of the panels it is clear that star-formation 
cannot excite the lines we see in NGC 1266. The top diagram 
shows that the gas is ionised either by an AGN, or shock processes. 
The middle and bottom panels show that the majority of the spaxels 
have line ratios which are consistent with low-ionisation nuclear 
emission region (LINER: iHeckmanl Il980h lik e activity. LINER 
are a controversial classification (see iHol 20081 for a review), with 
som e authors claiming that the ionisation comes from an AGN 



(e.g. iFerland & Netzer]|l983l:lHo. Filippenko & Sargentlll993l). fast 



shock s l lBurbidge. Gould & Pottaschl 1 19631 : iDopita & Sutherland! 
Il995h. or photoionization by ultraviolet (U V ) brigh t s tellar sources 
i di Serego Alighieri. Trinchieri & Brocatol Il990l ; iBinette et alj 
ll994l : ISarzi et alJl20ld : IO'Dea et alj|2010h . 



In addition to the BPT diagrams, it is also possible to plot the 
observed line ratios over grids that predict the line ratios of ionised 
gas lines in given scenarios. Both obscured and u nobscured AGN 
models from lGroves. Dopita & Sutherlan dl J2004I) fail to fully re- 
produce the observed line ratios. The closest AGN model requires 
a metallicity of 4Zq, an electron density of 1000 cm 3 , and a rel- 
atively constant radiation field. Even then this model cannot fully 
reproduce the observed line ratio distribution. We find that the only 
model that can repro duce the emission line ratios we observe is 
a shock model from lAUen et all d2008h . again with a super so- 
lar metallicity (« 2Zq). The stars in NGC 1266 have a sub-solar 
metallicity (McDermid et al., in prep), and thus if the gas really is 
of super-solar metallicity then it must have been significantly en- 
riched. 

Figure [14] shows our best fitting shock models, from the 
SAURON specific line ratio diagram [OIII]/H/i vs [NI]/Ha and the 
SAURON/GMOS combined [OIII]/H/3 vs [SII]/Ha and [OI]/H« 
diagrams. The emission observed requires shock velocities of up to 
~800km s -1 , and an average magnetic field of ~3 fiG cm 3 / 2 . The 
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shock velocities are similar to the outflow velocities, suggesting that 
the shocked emission is powered by the outflow. 

The shock velocities derived from these grids are highest at 
the centre, and towards the northern (redshifted) side of the out- 
flow (see Panel a of Figure[T5l>. This is true in all the data sets, even 
though the line ratio that drives the shock determination is different. 
The [SII]/Ha grid predicts higher shock velocities than the [NI]/H« 
grid, but the average shock velocity is consistent. The [OI]/Ha grid 
predicts lower shock velocities on average, but this is driven by the 
bins that fall within the Seyfert region of the associated BPT dia- 
gram (Figure[T3"l>. These spaxels are mostly located in the centre of 
the galaxy, where the AGN contribution to the ionisation is likely to 
be strongest. It is not surprising that suggestions of the embedded 
AGN are most clearly seen in the [OI]/Ha diagram, as the pres- 
ence of strong [OI]A6300 is usually indicative of a power-law ion- 
ising spectrum. The [OI] line, which is collisionally excited, will 
only occur in a zone which has a sufficiently high electron density 
and temperature to excite the upper level. With a stellar input spec- 
trum, these conditions only occur within the H + Stromgren sphere, 
where the 0° abundance is negligible. However, a gas ionised by a 
relatively flat power-law spectrum ha s an extended p artially ionised 
zone where the [OI] emission arises jPetersorJll997h . 

The magnetic field appears to increase with increasing distance 
from the centre of the galaxy, perhaps due to compression of mag- 
netic field lines as they enter the shocked gas (Panel b; Figure [T5V 
Shock velocities are similar in both the systemic and outflow com- 
ponents, again arguing that these components are linked, but the 
magnetic field may be higher, on average in the outflowing compo- 
nent. 

A2011 have shown that NGC 1266 has a deeply embedded 
AGN, but as we show here, the ionised gas emission in this LINER 
galaxy is dominated by shocks, likely caused by the outflow. A sim- 
ilar shocked, high metallicity wind has been postulated to explain 
the L INER like emission from starburst-superwinds (e.g. lRich et alj 
l2010h . The shocks in NGC 1266 exist at high velocities, similar to 
those found in the outflow, once again demonstrating the link be- 
tween the ionised gas and the material in the outflow. 

3.4 Physical Conditions 

Using the ionised gas lines we detect, it is possible to estimate the 
physical conditions within the ionised gas in each spaxel. From 
the ratio of the 6731 A and 6717 A lines in the [SII] doublet it is 
possible to estimate the luminosity-weighted electron density along 
each line of sight. The ratio of the [Nil] doublet at A6548,6584 A 
and the [Nil] line at A5755 A also provides a sensitive tempera- 
ture tracer, that depends very little upon the density of the gas. We 
calculated these line ratios for each spaxel where it is constrained 
by our GMOS data, and used the nebular package in IRAF to es- 
timate the electron density and temperature, based on the 5-le vel 
atom program developed bv lDe Robertis. Dufour & Huntl Jl987h . 

Figure [T6l shows maps of the derived electron density for the 
systemic and outflow component, and a histogram of the bin values 
in both components. In significant areas of the outflow component 
the [SII] ratio is saturated at the low density limit. As shown in 
panel b of Figurell6l however, in the highest velocity regions of the 
outflow where the line ratio is not saturated, the average density is 
higher than that in the systemic component. The required lines to 
calculate the electron temperature were detected in fewer spaxels, 
and so we do not reproduce these plots here. However a reasonably 
consistent picture emerges, with the higher density outflow having 
lower average temperature, and vice versa for the lower density sys- 
temic component. We hypothesize that the outflow may be higher 
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Figure IS. Shock velocities (top panel) and magnetic field strengths (bot- 
tom panel) for the systemic component derived from the SAURON only 
[OIII]/H/3 vs [NI]/Ha diagram by interpolating the shock grid presented in 
Figure [T4l The same trends are observed in the mixed GMOS/SAURON 
diagrams. 

average density due to a pile up of material in-front of the ionised 
gas shock, and/or an increase in the number of free electrons, due to 
ionisation from the shock. The disturbed structure in the systemic 
component appears to have a lower electron density than the rest of 
the gas, on average, supporting the hypothesis that this could be a 
different gas component. 

The northern half of the ionised gas systemic component has 
a lower electron density, on average, than the southern half. Some 
of the ionised gas densities in the northern regions are even at the 
low density limit. This could potentially offer an explanation for 
some of the puzzling features of the outflow. If the ISM is unevenly 
distributed on each side of the nucleus, then the same input energy 
from the central engine would dissipate faster in the denser material, 
producing stronger line emission in the denser southern regions, as 
observed. The lower density in the northern region would allow the 
same energy input to produce a faster shock, as observed. The mag- 
netic field in the northern region may also be lower, and this would 
decrease the expected synchrotron flux from a jet, as observed in 
Section [3~5l One might however expect that the outflow would thus 
be faster on the northern side of the nucleus, as the out-flowing ma- 
terial would encounter less resistance. This is not observed in the 
velocity fields presented in this paper. 

The observable tracers of the outflow that we currently posses 
(ionised and atomic gas in this paper, molecular gas in A201 1 ) how- 
ever have to be present at sufficient densities to be detectable, and 
more sensitive observations may be required to detect the high- 
est velocity parts of the redshifted outflow. Both evidence from 
A2011, and from the NaD absorption presented in this paper sug- 
gest that the redshifted northern part of the outflow is seen through 
the galaxy, and dust obscuration could potentially prevent us from 
detecting low level ionised gas emission. 
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Figure 16. Panel a shows maps of the electron density (in log units), derived from the [SII] doublet, within the systemic component (left) and outflow (right 
panel) of NGC 1266. Panel b shows a histogram of the electron density from the non-saturated bins, for the systemic component (solid line) and the outflow 
(dashed line). The median error in the electron density is shown in the legend. 



3.5 Driving mechanisms 

The physical mechanism by which an AGN removes gas from a 
galaxy is still widely debated. As discussed above, radiation pres- 
sure, broad-line winds, cosmic ray heating, and radio jets have all 
been mooted as potential methods for driving gas out of galaxies. 
In NGC 1266 the gas is clearly being accelerated and removed from 
the galaxy, and the cold gas phase survives reasonably far out into 
the outflow. A2011 argue that radiation driven winds are unable to 
provide enough mechanical energy to drive the outflow, and thus an 
AGN jet is a more likely driving mechanism. 1 .4 Ghz and 5 Ghz ra- 
dio co ntinuum maps are available for NGC 1266 dBaan & Klocknej 
l2006h . and are shown here in Figure [T7] The 1.4 Ghz image shows 
a central bright point source, and a jet-like extension to the south 
and north. 

A similar jet-like structure is visible in the 5 Ghz image, but 
in this higher resolution image the emission to the south and north 
of the nucleus is not connected to the point source. This lack of 
connection could be real, or because the VLA observations used 
here have resolved out the connecting structure. We could be seeing 
emission from where the jet hits the surrounding ISM (c.f. hotspots) 
or the jet itself may contain discrete blobs of synchrotron emitting 
plasma. 

The central source in this object appears to be driving the out- 
flow, despite being a low-luminosity AGN (having a total radio 
power at 5 GHz of < 10 22 W Hz -1 ; Alatalo et al., in prep). This 
highlights that low-luminosity AGN may also be important sources 
of feedback. The spectral index of the central point source (which 
is unlikely to be affected by itu-coverage issues) has a spectral in- 
dex typical of an low-luminosity AGN source (a[1.4-5Ghz]~-0.66; 
Alatalo et al., in prep). The spectral index when including the outer 
structures steepens (a[1.4 -5Ghzlsa-0.79), as is seen in ho tspots in 
local radio galaxies (e.g. Irloekstra. Barthel & Hej Il997t) . but as 
some emission could have been resolved out, we cannot be sure 
this value is robust. 

Figure [T8l shows the 1.4 Ghz radio continuum structure over- 
plotted on our GMOS maps of the ionised gas and neutral gas out- 
flow. The observed radio structure correlates well with the features 
observed in the ionised and atomic gas. The radio morphology is 



similar to that predicted by models of disrupted Fanaroff and Riley 
Class I (FRI) objects, as expected given the very dense environment 
around the jet, where it can be rapidly decelerated and disrup ted by 
gas entrainment (e.g. lLaing & Bridlell2002l : lLaing et~al]|201ll) . 

The 'kidney-bowl' shock structure observed in the ionised gas 
traces well the southern edge of the jet seen in radio continuum, 
as is observed in many radio galaxies (i.e. the "Alignment E ffect", 
IChambers. Milev & van Breugelll987l : lMcCarthv et al.ll987l) . This 
suggests that the shock suggested by the ionised gas line ratios oc- 
curs where the radio plasma impacts with the ionised gas. The blue 
shifted feature on the south-eastern side of the systemic component 
also seems to correlate well with the jet. This feature may well be 
associated with the jet, arguing that the systemic component may 
not represent a single structure. 

The radio continuum structure hooks around the fastest mov- 
ing neutral and ionised gas. As discussed above, the position angle 
of the outflow traced in NaD (and molecular gas; see Figure fTOb 
is slightly different than that traced by the ionised gas. Projection 
effects could possibly explain some of the difference between the 
position angles, but is unclear if this can completely remove the dis- 
crepancy. If this effect is real, then it could potentially arise because 
interaction with the ISM has deflec ted the jet, or the jet could be pre - 
cessing (as observed e.g. in Cen A: lHavnes. Cannon & Ekersll983l) . 
Of course it is also possible that the jet itself is powering the ionised 
gas outflow, but that another mechanism (e.g. winds or radiation 
pressure) could be driving the atomic and molec ular outflows, natu- 
rally leading to different position angles (see e.g lRupke & Veilleuxl 
l201ll) . From the speed of the outflow, and its measured extent we 
can estimate the lifetime of the phenomenon. Both the ionised gas 
(extending to ~ 1.5 Kpc and travelling at upto 900 km s _1 ) and the 
neutral gas (extending ~400 pc and travelling at an average 250 km 
s -1 ) give a similar estimate for the lifetime of the outflow (~1.6 
Myr). This would seem to suggest the outflows in ionised and neu- 
tral gas have a common cause. Radiation pressure driving the out- 
flow also seems unlikely on energetics grounds (see A2011). Deep 
high resolution observations of the radio continuum structure could 
help distinguish the cause of this effect. 
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Figure 17. 1 .4 Ghz and 5 Ghz radio continuum maps of the centre of NGC 1266 from lBaan & K lockner 1 200(J). The maps are shown for all emission above 
5tr (1.4 Ghz) and 3cr (5 GHz) with contours spaced by Icr. The beam sizes of sal'.'S and Rify/4, respectively, are indicated in the figures in the bottom left. 
NGC 1266 has an integrated flux density of 90.34 ± 2.83 mjy at 1.4 Ghz (Alatalo et al., in prep). 
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Figure 18. 1.4 Ghz radio emission contours (spaced every 8% from 90% to 8% of the peak flux) over plotted on Ha flux maps for the systemic and the outflow 
components (panel a & b), velocity fields for the systemic and the outflow components (panels d & e) and NaD neutral gas EW and kinematics (panels c & f ). 
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4 Conclusions 

In this paper we have been able to shed some light on the out- 
flow activity in NGC 1266. This unusual galaxy is relatively nearby, 
allowing us to investigate the process of AGN feedback in action. 
Using the SAURON and GMOS IFUs we detected strong ionised 
gas emission lines (Ha, H/3, [OIII], [01], [Nil], [SII] and Hel), as 
well as NaD absorption. We use these tracers to explore the struc- 
ture of the source, derive the ionised and atomic gas kinematics and 
investigate the gas excitation and physical conditions. 

NGC 1266 contains two ionised gas components along each 
line of sight, tracing the ongoing outflow and a component closer to 
the galaxy systemic, the origin of which is unclear. The gas appears 
to be being disturbed by a nascent AGN jet. 

We have presented further evidence the outflow in this object 
is truly multiphase, containing radio emitting plasma, ionised and 
atomic gas, as well as the molecular gas and X-ray emitting plasma 
(as detected in A201 1). With outflow velocities up to 900 km s , 
the outflow is very likely to be removing significant amounts of 
cold gas from the galaxy. The ionised gas morphology correlates 
well with the radio jets observed in 1.4 Ghz and 5 Ghz continuum 
emission, supporting the suggestion of A2011 that an AGN jet is 
the most likely driving mechanism for the ionised gas outflow. 

The line emission in NGC 1266 causes it to be classified as a 
LINER in optical diagrams. We show here that although NGC 1266 
undoubtably hosts an AGN (see A2011) the line emission in this 
object is extended, and is most consistent with excitation from fast 
shocks caused by the interaction of the radio jet with the ISM. These 
shocks have velocities of up to 800 km s~ , which match well with 
the observed velocity of the outflow. 

Using the observed NaD absorption we are able to set fur- 
ther constraints on the size and orientation of the outflow. We show 
that we are able to detect atomic gas entrained in the outflow out 
to a (deprojected) distance of ~400±50 pc, and that the outflow 
has an inclination (between the galaxy plane and the outflow) of 
53±8°. This suggests the outflow is misaligned from the stellar 
body. Furthermore we were able to provide an independent es- 
timate of the column density of neutral material in the outflow 
Nhi=(1.2±0.6)x 10 21 cm -2 . This estimate is consistent with that 
derived from HI absorption in A2011. The neutral and molecular 
outflows are well correlated, but appear to be outflowing along a 
slightly different axis to the ionised gas. The cause of this affect is 
unclear. 

NGC 1266 is a highly complex object, and it is clear that fur- 
ther observations will be required to fully understand it. Observa- 
tions of single emission lines (either with an IFU or with a Fabry- 
Perot instrument) will be important to overcome problems with line 
blending inherent in this high velocity dispersion source. Higher 
spatial resolution would also be advantageous to obtain clear di- 
agnostics of AGN activity, and better understand the shock struc- 
ture. Sensitive interferometric radio observations at high angular 
resolution would also enable us to understand the morphology and 
orientation of the nascent radio jet. The Atacama Large Millime- 
ter/submillimeter Array (ALMA) will allow us to study the molecu- 
lar component of the outflow in greater detail, and further determine 
the driving mechanism. For instance if molecular shock tracers are 
found predominantly in the outflow then a kinetic driving mecha- 
nism would be favoured, while if photon dominated region tracing 
species were detected this would argue for a radiation driven com- 
ponent to the outflow. 

This galaxy is one of the few currently known in which we can 
witness ongoing active feedback, where a central AGN is disrupt- 



ing its star-forming reservoir. It is clear that understanding the pro- 
cesses removing the ISM will have widespread applications to both 
theoretical and observation attempts to understand AGN feedback, 
its affect on the ISM, and role in building up the red-sequence. 
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